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Abstract 

We discuss a theory of the ultrafast non-linear optical response of excitons in the 
Quantum Hall system. Our theory focusses on the role of the low energy collective 
electronic excitations of the cold, strongly correlated, two-dimensional electron gas 
(2DEG) present in the ground state. It takes into account ground state electron 
correlations and Pauli exchange and interaction effects between the photoexcited 
excitons and the collective excitations. Our formulation addresses both the initial 
coherent regime, where the dynamics is determined by exciton and 2DEG polariza- 
tions, and the subsequent incoherent regime, dominated by population dynamics. 
We describe non-Markovian memory, dephasing, and correlation effects and a non- 
linear exciton hybridization due to the non-instantaneous interactions. We iden- 
tify the signature of the coherent inter-Landau level magnetoplasmon (MP) in the 
temporal profile of the three-pulse time-integrated four-wave-mixing signal in the 
initial coherent regime. 
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1 Introduction 



In low dimensional geometries, collective many-electron effects can lead to 
new quantum phases that display novel transport and optical properties. One 
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such example is the cold two-dimensional electron gas (2DEG) in modulation- 
doped quantum wells (MDQW), which provides an ideal "laboratory" for 
studying and manipulating charge and spin excitations on the nanometer and 
femtosecond scale [1,2,3]. In the presence of a strong magnetic field point- 
ing along the MDQW growth direction, the 2DEG is a strongly correlated 
quantum liquid whose ground state cannot be described by a non-interacting 
Hartree-Fock state, unlike the ground state of undoped semiconductors [1,2]. 
The magnetic field creates discrete Landau levels (LL), which in the ground 
state are partially filled by electrons introduced by the doping. These LLs 
are highly degenerate, and the states of each LL are strongly coupled by the 
Coulomb interaction. The latter coupling leads to collective charge and spin 
excitation modes, such as the intra-LL and inter-LL magnetoplasmons (MP), 
with dispersion governed by strong correlation effects that lead to a magne- 
toroton minimum [1,2]. The LL degeneracy increases with magnetic field, and 
above a threshold value the ground state electrons only occupy the lowest LL 
(LLO) states; all the higher LLs (LL1, LL2, • • • ) are then empty in the ground 
state. The percentage of occupied states gives the LL filling factor, v. The 
ground state, well described by the Jastrow-like Laughlin wavefunction [1,2], 
depends on this filling factor. Strong exchange Coulomb interactions stabilize 
a ferromagnetic ground state with polarized electron spins for filling factors 
v — 1/m, where m is an integer [2,4], or for integer v. For fractional values of 
v, the correlated ground state leads to the fractional QHE [1,2,3]. 

Recent time-resolved four-wave-mixing experiments shed new light into the 
dynamics of this strongly correlated system, and opened a new field of non- 
equilibrium Quantum Hall physics [5,6,7,8,9,10,11,12,13,14]. The interband 
optical absorption spectrum of the 2DEG is dominated by LL exciton peaks 
[5,8]. In the QHE system, the 2DEG correlations and collective electronic ex- 
citations dominate the dephasing of the photoexcited excitons (X) for low 
photoexcitation intensities [5,6,7,8,9,10,11,12,13,14]. Of particular interest is 
the dynamics of the QHE system during time scales comparable to the char- 
acteristic time it takes the cold 2DEG "bath" system to react to the introduc- 
tion of photoexcited Xs. The dephasing times of the two lowest LL excitons 
in the QHE system range from a few picoseconds (LLO) to a few hundreds of 
femtoseconds (LL1), while the 2DEG responds to the X-2DEG interactions 
within a time interval comparable to the period of its low energy excitations. 
The period of the lowest inter-LL magnetoplasmon (MP) collective modes 
[1,2], discussed below, is Tmp = IttH/Qm, where Qm ~ 15 — 20 meV is the 
MP excitation energy. Therefore, Tmp is of the order of a few hundreds of 
femtoseconds, longer than the duration of the ~100fs optical pulses used to 
probe this system and tunable by changing the magnetic field. Strong quan- 
tum kinetic effects in the ultrafast non-linear optical dynamics are expected 
in this regime, explored for the first time in [5,6,7,8,9,10,11,12,13,14]. In order 
to study time scales comparable to the duration of the X-2DEG scattering 
events, well-established pictures such as the semiclassical Boltzmann picture 
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of dephasing and relaxation (which assumes instantaneous X-2DEG scatter- 
ing) and the Markov approximation need to be revisited [15,16,17,18,19,20]. 

Transient wave-mixing spectroscopy is well suited for studying strong cor- 
relation and quantum coherent effects during time scales shorter than the 
dephasing times. Indeed, fluctuations beyond the Random Phase Approxi- 
mation (RPA) can generate a FWM signal with a distinct time dependence 
[15,16]. For example, in undoped semiconductors, exciton-exciton (X-X) in- 
teractions dominate the two-pulse FWM signal during negative time delays, 
where the Pauli blocking effects vanish [15,16]. The time-dependent Hartree- 
Fock (HF) treatment of such interactions [21] predicts an asymmetric FWM 
temporal profile, with a negative time delay FWM signal decaying twice as fast 
as the positive time delay signal [15,16,22]. The observation of strong devia- 
tions from the asymmetric HF temporal profile is the signature of strong X-X 
correlations [15,16]. In this paper, our main goal is to describe theoretically 
the non-instantaneous many-body effects in the ultrafast non-linear optical 
response of the QHE system that are due to the correlated 2DEG. We discuss 
ways in which ultrafast non-linear optical spectroscopy can be used to probe 
the coherent MP dynamics, important for understanding and manipulating 
the coherent transport and optical properties of the 2DEG but not accessible 
so far. 

The two-pulse wave mixing spectroscopic studies of electronic dephasing in 
the QHE regime reported in [5,6,7,8,9,10,11,12,13] focussed on the dynamics 
of the excitonic peaks for filling factors v < 2, where only LLO states are pop- 
ulated in the ground state. The time and frequency dependence of the FWM 
spectra revealed new dynamical features that could not be explained within 
the RPA. The most important observations were a strong optically-induced 
time-dependent hybridization of the LLO and LL1 exciton peaks governed by 
the inter-LL MP dynamics, a non-Markovian dephasing of the LL1 exciton 
peak when v < 2, and strong oscillations as function of time delay of quantum 
kinetic origin. The exciton hybridization manifests itself in a strong enhance- 
ment of the LLO FWM signal even in the case of minor photoexcitation of the 
LLO optical transitions. This non-resonant FWM signal displayed a symmetric 
time delay profile, unlike for the LL1 signal and the predictions of the RPA. 
Additionally, a strong dephasing of the LL1 resonance and an asymmetric 
LL1 linear absorption lineshape was observed when v < 2, unlike for the LLO 
resonance. Finally, for equal photoexcitation of the LLO and LL1 transitions, 
the FWM signal at the LL1 energy was very minimal, while the LLO FWM 
resonance showed strong oscillations as function of time delay, despite the fact 
that a single LLO peak dominated the FWM spectral profile. In contrast to the 
MDQW, a comparable undoped quantum well did not display any of the above 
behavior: its FWM spectrum was consistent, to first approximation, with the 
predictions of a three level system consisting of the ground state, LLO, and LL1 
levels. Importantly, the strength of the new MDQW features diminished as the 
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density of photoexcited carriers approached the density of the ground state 
2DEG. In this case, the Coulomb scattering among the photoexcited carriers 
dominates over the X-2DEG interactions and can even change the dephasing 
times. More recent experiments by the Chemla group [14] showed that the 
initial coherent regime, which lasts for a few picoseconds, is followed by an in- 
coherent regime that, for low photoexcitation intensity, builds up slowly over 
a time interval of several picoseconds. 

To interpret experiments such as the above, one should note some important 
differences between the MDQW and QHE system and the undoped QWs stud- 
ied so far. In addition to the ground state electron correlations in the QHE 
system, important differences in the optical properties between the doped and 
undoped QWs arise from the different nature of the low energy electronic 
excitation spectrum in the two systems. In undoped semiconductors, the low- 
est electronic excitations of the ground state electrons are the high energy 
interband e-h pair excitations, which can adjust almost instantaneously to 
the dynamics of the photoexcited carriers. The latter then behave as quasi- 
particles with mutual interactions, while the ground state can be considered 
as rigid: the many-body nature of the system only affects the different param- 
eters associated with the band structure and the dielectric screening. 

In the case of the 2DEG in the QHE regime one should distinguish between 
the excitations of two subsystems: (i) the QW interband excitations (with the 
2DEG at rest), which consist of 1 e-h, 2 e-h, ■ ■ ■ pairs created in the LLs, and 
(ii) the 2DEG excitations (with unexcited QW and full valence band), e.g. 
the 1-MP, 2-MP, • • • and incoherent pair excitation states. The ensemble of 
states that determine the non-linear optical spectra can be thought of as con- 
sisting of £ e-h pairs and n 2DEG excitations. One can then draw an analogy 
between the X-MP effects studied here and the X-phonon effects in the ultra- 
fast response of undoped semiconductors that were extensively studied in the 
past [17,18,20,23,24,25,26,27]. In both cases, the ultrafast non-linear optical 
dynamics is governed by the interactions between photoexcited Xs and a col- 
lective excitation. However, there are some important differences. In the QHE 
system, the collective excitations are electronic in nature, and are described 
by similar electronic operators as the photoexcited carriers. They are therefore 
subject to Pauli correlations with the photoexcited Xs, while the ground state 
electrons are strongly correlated and sustain low-energy excitations. On the 
other hand, in the undoped system, the electronic operators commute with the 
collective excitation (phonon) operators, while the ground state correlations 
can be neglected. 

The role of the Coulomb and exciton-phonon correlations in undoped semi- 
conductors have been studied in the past using different theoretical approaches 
[18,20,23,24,28,29,30,31,32,33,34,35]. A widely used approach is the "dynamics- 
controlled truncation scheme" (DCTS) [19,23,24], which truncates the hierar- 
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chy of density matrices generated by the interactions based on the fact that, in 
the undoped system, all Coulomb interactions occur between photoexcited e- 
h pairs and are thus dynamically generated by the optical excitation, treated 
with an expansion in terms of the optical field. The DCTS also assumes the ab- 
sence of free carriers in the ground state [23], a condition that does not hold in 
the MDQWs. The almost unexplored dynamics of strongly correlated systems 
whose ground state electrons interact unadiabatically with the photo-excited 
e-h pairs raises very fundamental questions. In MDQWs, the direct exciton- 
exciton interactions are screened, and the non-linear response is mainly deter- 
mined by the Fermi sea excitations [34,36,37,38,39,40,41]. In the QHE regime, 
the presence of collective low energy electronic excitations and the resulting 
non-Markovian dynamics and memory effects, as well as the strongly corre- 
lated ground state, raise formidable theoretical difficulties for describing the 
non-linear optical dynamics. For example, standard diagrammatic expansions 
and DCTS factorizations that assume a Hartree-Fock reference state break 
down. In addition, the theory must address the quantum effects due to the 
Pauli correlations between the collective electronic excitations and the pho- 
toexcited carriers. 

In the first part of this paper we discuss a theoretical formulation that de- 
scribes the ultrafast third-order non-linear optical response of the QHE system 
at zero temperature and addresses both the coherent and incoherent regimes. 
To our knowledge, this is the first theory that addresses this non-equilibrium 
many-body problem. Our approach is based on the projection of the exciton 
states and the separation of the uncorrelated contributions to the third-order 
non-linear optical response from the contributions due to correlations among 
the interband and intraband elementary excitations. We also introduce a basis 
of correlated states, generated by using the Lanczos iterative approach [42], 
which is useful for calculating the 2DEG correlations [1,2] as well as for de- 
riving a simple average polarization model [16,30,35] for treating the ultrafast 
dynamics. We neglect phonons since the optical phonon frequency is much 
larger than the inter-LL splitting. In the second part of the paper, we focus 
on the coherent regime and discuss the manifestations of the MP coherence in 
the three-pulse time-integrated FWM signal, while neglecting the incoherent 
effects. This calculation is based on an average polarization model derived 
from the general theory by projecting out the first few correlated Lanczos 
states. 

Although our theory can be used to treat the general case, of particular inter- 
est here is the spin-f polarized LLO 2DEG, realized for filling factors v — l/m, 
where m is an integer, or integer v. We consider the case of photoexcitation 
with right-circularly polarized optical pulses, which excite spin-j electrons. 
For simplicity we only consider the LLO and LL1 electron and hole states. 
We show that the MP effects can be separated from the corresponding Pauli 
blocking and exciton-exciton interaction signals based on the temporal pro- 
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file of the FWM signal. Throughout the paper we point out the analogies 
between the theoretical formulation of the X-MP effects presented here and 
the DCTS [23,24] and correlation expansion [20] approaches for treating X- 
phonon effects in undoped semiconductors. Similar to the DCTS [24], we use 
an expansion in terms of the optical field in order to eliminate the number 
of independent dynamical variables that need to be considered. The X-2DEG 
correlations do not allow the complete factorization of the intraband density 
matrix into products of interband coherences. Here we separate out the cor- 
related contributions, which lead to the incoherent effects, without assuming 
a Hartree-Fock ground state, as in the case of the DCTS. In undoped semi- 
conductors, our approach reduces to the DCTS if phonons are included. 

In Section 2 we set up the general problem and in Section 3 we discuss the 
nature of the exciton and magnetoplasmon states that determine the optical 
spectra. We also present some useful relations for describing the Pauli ex- 
change quantum effects and the interaction effects between Xs and MPs. In 
Section 4 we present the equations of motion for the non-linear polarizations 
and photoexcited carrier populations and identify the contributions due to 
the many-body correlations. In Section 5 we present a decoupling scheme for 
treating the interaction effects, which is motivated by a decomposition of the 
photoexcited many-body states that separates out the uncorrelated and exci- 
tonic contributions from the correlated and incoherent contributions. We use 
this approach to devise a factorization scheme and identify the intraband and 
interband correlated contributions to the density matrix. We also discuss the 
linear absorption spectra. In Section 6 we briefly discuss the coherent X-X 
interaction and scattering effects and derive an average polarization model 
[16,30,35] for treating such correlation effects by using a basis of interacting 
X-X Lanczos states [42]. In Section 7 we derive from the general equations a 
generalized average polarization model that treats the effects of the magneto- 
plasmon dynamics. Finally, in Section 8, we use the above model to calculate 
the three-pulse time-integrated FWM signal and identify the signatures of 
the MP dynamics in the coherent regime by assuming that incoherent effects 
can be neglected for such short time scales. We end with the conclusions. 



2 Hamiltonian 

To describe the optical response of the QHE system, we adopt the standard 
two-band Hamiltonian of interacting electrons and holes coupled by an optical 
field E(t) and subject to a magnetic field, treated within the Landau gauge 
A = (0,Bx, 0), that splits the conduction and valence bands into discrete 
electron (e) and hole (h) LLs [21,29]: 

H tot (t) = H- fiE(t)Xi - fiE*(t)X, (1) 
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where X is the interband optical transition operator, discussed further in 
the following section, \x is the interband transition matrix element and N = 
L 2 /2rrl 2 is the LL degeneracy, where / = (hc/eB) 1 ^ 2 is the magnetic length 
(Larmor radius) and L the system size [35,43]. H is the "bare" semiconductor 
many-body Hamiltonian (h — 1 from now on) [21,29,35,7], 

H = J2\E g + n c c (n + 1/2)] e{ na e kna + ]T QF e (n + 1/2) h{ na h kna + H rnt , (2) 

kncr kncr 

where E g is the bandgap and VL l c = eB/rrii, % = c,v, are the electron and 
hole cyclotron energies that determine the Landau level spacings, which are 
inversely proportional to the carrier masses and closely spaced in the valence 
band. In the above Hamiltonian, e\ na is the creation operator of the spin a 
LLn conduction band electron state 

e iky 

^fcn(r) = -Jj$n{x ~ X k ) (3) 

characterized by the momentum k = x k /l 2 , where x k is the x coordinate of 
the cyclotron orbit center. In the above equation, <p n is the n-th eigenfunction 
of the ID harmonic oscillator, with energy equal to the cyclotron energy. 
Similarly, h kncr creates the LLn valence band hole state ■0, which in the ideal 
two-dimensional system is related to the conduction electron wavefunction by 
^kn = ip-kn- We only consider heavy hole states, with angular momentum 
J = 3/2 and m = 3/2 (a =|) or m = —3/2 (a =|). The above relation 
between the conduction electron and valence hole wavefunctions in the ideal 
two-dimensional system leads to an electron-hole symmetry that strongly 
affects the non-linear optical properties [35,44]. In the realistic system, this 
symmetry is lifted due to the lateral confinement, the different band offsets 
and confinement between the electrons and the hole, the valence band mixing, 
etc. 

The Hamiltonian H int describes the e-e, e-h, and h-h Coulomb interactions, 



Hi n t „ 



v ee < st p f e e +v hh ft* fit h h 



- v eh p f e h - v he P h e 



, (4) 



where a = (k,n,a). In the ideal two-dimensional system, the Coulomb inter- 
action matrix elements v l J ia2 a30l4 (with i,j = e, h) are given by 

V ct 1 a 2 ,a 3 a 4 ~ J J^^ V q^a 1 a 2 ( ( i) ■Ea 3 a 4 (~ ( i) ^ (5) 

where v q = 2ne 2 /q is the 2D Coulomb potential, and [43] 
^ e 1Q2 (q) = VnMe^+^Sto^S^ , i£ loa (q) = Ff Q2! _ m (q), (6) 
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where —a = (—k,n,cr) and 

L:-"(^f)e-«" 2,i , (7) 

for m > n and ip mn (<l) = ^Pnm('~ c l) f° r m < n i where L™~ n is the generalized 
Laguerre polynomial. From now on we measure all energies with respect to 
the ground state energy and thus have that H\G) = 0. 



<Pr. 



n\ 



-% + iq x )l 
V2 



3 2DEG Magnetoexcitons and Magnetoplasmons 



In this section we briefly discuss the magnetoexciton (X) and magnetoplasmon 
(MP) excitations that govern the ultrafast non-linear optical dynamics. From 
now on we restrict to the case of photoexcitation with right-circularly polar- 
ized light, which excites spin- J, e-h pairs. We start with the dipole transition 
operator X\ which can be expanded in terms of the exciton creation operators 
JtJ that create the allowed optical transitions for given LLs. We introduce the 
creation operators of LLm — > LLn magnetoexcitons with total momentum q, 

Vt — 1 X^ p ikq x l 2 -t U (o\ 

A q«m ^/jyZ^ C °k+q y /2,n,l u -k+q y /2,m,{- \°) 

In the absence of disorder, momentum is conserved and only q = excitons 
are photoexcited directly. Furthermore, in the ideal system, the only allowed 
optical transitions correspond to m — n. We then have that 

X^ = VNJ2XI, Xl = Xl n . (9) 

n 

We note however that the disorder present in the MDQW system can relax 
the momentum conservation condition [45,46,47,48], thus mixing X states with 
different momenta, while the valence band mixing couples the n ^ m valence 
hole states and magnetoexcitons [45,46,49]. For this reason, in the formulation 
of the non-linear optical response we first consider all allowed magnetoexciton 
transitions, X^ = VNJ2iX}, before presenting results specific to the ideal 2D 
system. 

The states \Xt) = X\\G) are the magnetoexciton states in the 2DEG system. 
The difference from undoped semiconductors is that here the exciton operators 
Xj act on the strongly correlated state \G), which is the ground eigenstate of 
the many-body Hamiltonian H that describes the correlated 2DEG at rest. 
Therefore, the states \Xj) are strongly correlated, similar to the magnetoplas- 
mon states discussed below. The following orthogonality relation holds: 

{X v \Xi) = {l-Vi)5 w , (10) 



8 



where, as discussed below, vi is related to the LL filling factor. For the magnetic 
fields of interest, only the electron LLO is partially filled in the ground state 
with the 2DEG at rest. On the other hand, all the hole LL states are empty 
(full valence band), similar to the undoped system. 

We now turn to the magnetoplasmon modes, which dominate over quasi- 
electron - quasi-hole pair excitations for momenta q < l/l. A LLm — >LLn 
MP may be thought of as an e-h pair, or exciton, formed by an electron in 
LLn and a hole in the LLm 2DEG. The creation operator of this MP is given to 
first approximation by the LLm — >LLn contribution to the collective density 
operator [1,2,50,51]: 

Pqnma = X! e ^ ^k+q y /2,n,a &k-q y /2,m,a- (H) 

The analogy between the above MP and magnetoexciton creation operators 
is clear. It is convenient to also introduce a similar collective operator for the 
hole states [44], 

Pqnma = X! C ^ ^-k+q y /2,n,a h-k-q y /2,m,a, (12) 

and note the relation 



N k 



Pqnma ~ P-qmna ? 2 — e, /i (13) 

that relates the creation and annihilation operators. Here we focus on photoex- 
citation of the LLO and LL1 optical transitions only, which are dynamically 
coupled by the LLO -> LL1 inter-LL MPs [5,7,8,9,10,11]. Theese MPs are the 
lowest-energy neutral charge excitations of the v — 1 QHE ferromagnet, where 
the intra-LL charge excitations are suppressed since all spin-f LLO states are 
occupied in the ground state [4]. 

Similar to Feynmann's theory of the collective charge excitation spectrum of 
liquid helium [1,2,50,52], a good variational approximation of the MP eigen- 
states of the Hamiltonian H is given by the state (single-mode approximation) 
[1,2] 

|M q ) = £ C nmCT (q)0 nm (q)^ nma |G), (14) 

anm 

where C nma (q) are variational parameters. The mixing of the different LL 
states, which is due to the interactions, is suppressed in the strong magnetic 
field limit [2] by a factor oc B^ 1 / 2 since the characteristic 2DEG Coulomb 
interaction energy e 2 /l oc \[B is smaller than the energy separation between 
the electron LLs, Vt c c oc B. Even though in the realistic system e 2 /l ~ 
calculations have shown that the LL mixing does not change qualitatively the 
MP and fractional QHE properties [43,50]. 
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We now turn to the Pauli exchange effects between the Xs. As already known 
from undoped semiconductors, the latter are described by the deviation of the 
commutator of the X operators from bosonic behavior due to the underlying 
Fermi statistics. By using the second quantization expressions of X and p and 
Eq. (10) we obtain that 

XqnmiX^, = (X qnm |X' t ) — — = A/31 qn / n j_ — — =/3^ qn , m |, (15) 



where Ap qnm | = p qnm j — (G\p^ nm ^\G) and, for the strongly correlated 2DEG 
ground state, the latter expectation value vanishes for nonzero momentum 
(and for n, m ^ in the strong magnetic field limit). More generally, we have 
that 

[X i ,Xj] = (J w (l-i/ i )-Ai> w (16) 

where the first term on the rhs describes the ground state contribution. In the 
case of zero-momentum excitons, 

[X n ,Xl} = 5 nn ,(l-v n -Av n ), (17) 

where 

^ = ^T,(G\el ni e kni \G) (18) 

k 

is the ground state filling factor of the LLn spin-j 2DEG electron states and 
1 

Az> ™ = Jj S (p-knfi-kni + e\ nl e k ni - {G\e\ ni e kni \G)) (19) 

k 

describes the change in the LLn filling factor due to the photoexcited electron 
and hole populations. 

Xs and MPs are made of electrons, and therefore the X-MP Pauli exchange 
effects must also be considered. Similar to the case of Xs, these are described 
by the commutators 

* t 1 "t h "t 1 "t 



obtained by using the above second quantization expressions for p and X . In 
the case of a spin-polarized ground state 2DEG of spin-f electrons, the MP 
and X operators commute since right-circularly polarized light creates spin-j 
electrons. 

In addition to the Pauli exchange effects, the optical properties are strongly 
affected by the interactions between the photoexcited excitons Xi and the 



10 



2DEG carriers. We can describe such X-2DEG interactions, which scatter the 
X into X+MP final states, by considering the action of the Hamiltonian H on 

\*i) [7]: 

H\Xi) = Q t \X t ) - (1 - Vi) ]T V Vi \X v ) + (21) 

The above equation defines the state \Yj) by the requirement that it is or- 
thogonal to all exciton states, (Xj\Yj) = 0, and therefore describes an excited 
2DEG configuration. Such configurations are denoted as 2DEG* from now on. 
The above orthogonality requirement, as well as the orthogonality among the 
exciton states, gives 

(22) 



the Xi energy, and 



v - {Xi\H\Xj,) 



the static Coulomb-induced coupling of the different LL Xs. Based on the 
above we introduce the operator 

Yi = [Xi, H] - n t X t + (1 - Vi ) ]T V W X V , (24) 

that describes the interactions between Xi and all the other carriers, X's, MPs, 
or ground state 2DEG. 

An explicit expression for the operator Y n was obtained in [7] by calculating 
the commutator [X n , H] using the Hamiltonian ifj nt discussed in section 2. 
By only retaining contributions from the photoexcited LLs (LLO and LL1) we 
obtain that 

[X n , H] = [E g + (n + 1/2) (Q c c + Q v c )] X n 

- (1 - V n ) X n> / -JTZ-yi V Q \4>nn'{q)\ 2 + (S n i ~ 5 n0 )Y int (25) 

where the first two terms give X energies and Coulomb-induced couplings 
similar to the undoped system [35,53], while the interaction contributions are 
described by the operator 

Yint = ~ ,} 7= v qPi [0io(-q)^ q oi - 0oi(-q)A q iol , (26) 

where we defined for simplicity 

Pq = drum' (<j) [p\mm>c ~ Pqm'ma) ■ ( 27 ) 



mm' a 
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The operator Y can be obtained from Eq. (24) by subtracting from the above 
expression for Yi nt the contributions to the X energies and couplings, Eqs. (22) 
and (23), which are oc (X n \Y^ nt \G) . If we restrict to the first two LLs we have 
the simple property 

Y = -% = Y. (28) 



In the undoped system, we have that Y^ nt \G) = and Y = Y int describes the 
XX interaction effects. In the doped system, Y int renormalizes the X energies 
and couplings due to the X interactions with the ground state 2DEG. In addi- 
tion to X-X interactions, Yi nt describes X+MP scattering effects. This can be 
seen from Eqs. (26) and (27) by recalling that the operators /5q mm / CT create and 
annihilate the MPs. To make the analogy to the case of phonons in undoped 
semiconductors [20,24], (Y) describes the contribution of the MP-assisted in- 
terband density matrices and the X-X interactions. These two contributions 
can be distinguished in the case of a spin-polarized 2DEG, where the MPs 
are excitations of the spin-f electrons that populate the ground state, while 
all a =1 carriers are induced by the right-circularly polarized optical pulses. 
In this case, the a =j contribution to Eq. (27) describes the MP interaction 
effects, while the a —[ term describes the X-X interactions. The results of 
[20,24] for the polarization equation of motion are reproduced if we add to 
the Hamiltonian H the electron-phonon interaction and use Eq. (24). The 
difference here is that both the X-X interactions and the interactions between 
the photoexcited carriers and the collective 2DEG excitations are described 
by the same electronic Hamiltonian Eq. (4). 

As demonstrated by Eq. (26), the state Y^\G) is a linear combination of {1- 
MP + 1-LLO-e + 1-LLl-h} and {1-MP + 1-LLl-e + 1-LLO-h} four-particle 
excitations, into which both the LLO and the LL1 excitons can scatter by 
interacting with the 2DEG. In the case of Xi, the LL1 electron can scatter to 
LLO by emitting a LLO — > LL1 MP. Since this MP energy is close to the e-LLO 
— > e-LLl energy spacing, the above scattering process is almost resonant. It 
therefore provides an efficient decay channel of the LL1 exciton to a {1-MP 
+ 1-LLO-e + 1-LLl-h} four-particle excitation of the ground state \G). All 
other allowed scattering processes are non- resonant. The Xi hole can scatter 
to LLO by emitting a MP, which leads to a {1-MP + 1-LLl-e + 1-LLO-h} four- 
particle excitation. The latter state however has energy that is significantly 
higher, by an amount of the order of ~ Vt c c + Vt v c) from that of the initial Xi 
state. In the case of X , the LLO electron can scatter to LL1 by emitting a MP, 
so that X — > {1-MP + 1-LLl-e + 1-LLO-h}, or the LL1 hole can scatter to 
LLO, in which case Xo — > {1-MP + 1-LLO-e + 1-LLl-h}. \Yq) is thus a linear 
combination of the same final states as |Y"i), also seen from Eq. (28). However, 
in this case the energy of all final states is significantly higher than that of 
the initial state \X ). Therefore, the decay of the LLO exciton is suppressed 
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as compared to that of the LL1 (or higher) exciton. As discussed below, this 
difference in the dephasing of the two LL excitons already plays an important 
role in the linear absorption spectra, and has even more profound effects on 
the non-linear optical spectra. 



4 Ultrafast Non-linear Optical Response 

In this section we obtain the equations of motion for the interband polariza- 
tions and the intraband populations and coherences that determine the non- 
linear optical response and identify the contributions due to the many-body 
interactions. 

Within the dipole approximation, the optical spectra are determined by the 
polarization of the photo-excited system, 

P (t) = -j=Y,m, * = (29) 

which may be obtained from the equations of motion for the Xi polarizations 
Pi. The time evolution of any operator O is determined by the full many-body 
Hamiltonian H tot (t), which here includes the 2DEG degrees of freedom that 
are responsible for the dephasing: 

id t {6) = ([d,H])- d(t) ]T<[6, x}\) - d*(t) £<[o, Jti]). (30) 

3 3 

where d(t) = f/,E(t)\/~N is the Rabi energy. Substituting O = Xi in the above 
equation and using the property [Xi, X^] = and Eq. (24) for the commutator 
[Xi, H] we obtain the familiar polarization equation of motion: 

id t Pi(t) - QiPiit) + (1 - Vi) Vn,Pi,{t) = -d(t)[l -Vi-Y. A ^'] + ^i), 

(31) 

where Az/^ = {AOw). The lhs of the above equation describes the static 
exciton energies and Coulomb-induced LL couplings. Two sources of non- 
linearity are described by the terms on the rhs and are the subject of this paper. 
The first term describes the Pauli blocking effects (PSF), which are determined 
by the ground state spin- J, electron populations, with filling factor and by 
the photocarrier populations and intra-band Raman coherences in the spin-j 
electron system, described by (Az>^), Eqs. (16) and (19). The second term on 
the rhs of Eq. (31), (Xi), comes from the interactions between Xi and the rest 
of the carriers in the system: X-X, X-MP, and X-2DEG interactions. 

The equation of motion for Aiv may be obtained by substituting O = [Xi, X\] 
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in Eq. (30) and recalling Eq.(16) that expresses Ai>ai in terms of the exciton 
operators. Calculating the commutator [H, Avai] by using the property 

[A, [B, C]] + [C, [A, B]] + [B, [C, A]] = 0, (32) 

which holds for any operators A,B,C, and Eq. (24) for the commutators 
[Xi, H], we obtain the equation of motion 

id t Av iif = - - ^) Av iV + {[Xi,Y}]) - ([Yi,Xj,}) 
+ (1 - v v ) Vji'Avij - (1 - Vi) VijAvji, 

- d(t)J2( [Av w , X}} ) - d* {t ) £ ( [Avu> , XA ) (33) 

j j 

where Ti is the relaxation time. In the ideal system, Av^i = Av^ai where 
Az/j, Eq. (19), is the photo-induced change in the LL filling factor. 

A simple model can be derived if we restrict to the first two LL states. Using 
Eqs. (28), (24), (32), and Eq. (17) for Av n we first obtain the commutator 

[A0 n , H] = (5 nl - 5 n0 ) (Ml - M n ) , (34) 

where we introduced the operator 

M n = [Y,Xl], (35) 

discussed below. Substituting the above into Eq. (30) and using the commu- 
tator 

[Av n ,XU = pn m Xl (36) 

obtained by using the second quantization expressions for Au n and X m , we 
obtain the equation of motion of the photo-induced change in the LLn filling 
factor, n — 0, 1: 



%d t Av n = - [d*(t)P n - d{t)P*\ + (<J n>1 - M [(M n y - (M n ) 



J 1 

(37) 



We note here that the above equation is consistent with the conservation of 
the total number of photoexcited carriers, Av + Ai/i. To see this, we first note 
the following property in the case of two LLs 

M n + Ml = V n , n [(1 - z/ )Az>! - (1 - vJAvo] , (38) 

where n' ^ n and n,n' = 0, 1, which can be obtained after some algebra by 
substituting in the definition Eq. (35) of M n the operator Y = Y\ = —Y 0: Eq. 
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(24), and using the commutator relation Eq. (32). Noting that Vio = Voi, we 
obtain from the above equation that M\ — M\ = Mq — M and 

id t (Ai/i + Au ) = ~ (Ai/i + Az/ ) + A [tT(*)(P + Pi) ~ d(t)(Po + Pi)] ■ 

(39) 

The above equation demonstrates the conservation of the total photoexcited 
carrier density when restricting to photoexcitation of the first two LLs. The 
relaxation time T\ is due to the scattering of the LLO and LL1 carriers to 
higher LLs as well as their radiative recombination, while the last term on the 
rhs describes the photoexcitation process. 

The intraband density matrices ([Y, X%\) that enter on the rhs of Eq.(37) 
describe a redistribution of the photoexcited carrier populations between the 
two LLs that is assisted by the MP and the interactions. Analogous phonon- 
assisted effects in the case of undoped semiconductors are discussed e.g. in [20] 
and [24]. The corresponding physical processes become clear by calculating the 
above commutator using Eqs. (26) and (8). Restricting for simplicity to the 
first two LLs we obtain that 

0io(q)(G|X qO1 Xt|G)-0 O i(q)(G|X qlo Xt|G)" 

+ (Snl-S n0 ) ^2 J2 [01o(q)A_ qO l-0Ol(q)A_ q io] [<j>10 (q) A_ q0 l - 001 (q) A_ q io 

+ [ 5 «o (0oi(q)p qOU -0io(q)p5oi|)-^i ( < Mq)p q ioi-0oi(q)Pqicn)_ 

(40) 

The first term on the rhs describes the photoexcitation of coherent MPs and 
is analogous to the coherent phonon contribution in undoped semiconductors 
[20]. Similar to the latter case, it vanishes in the ideal system, but is known 
to contribute in the realistic quantum Hall system due to disorder, inhomo- 
geneities, and valence band mixing [45,46,47,48]. Similar to the undoped sys- 
tem [24], the second term describes a contribution due to X populations and 
inter-LL exciton coherences, described by the density matrices (A" qnm A" qn / m /). 
The latter exciton coherences and populations come from the photoexcited 
carriers, while in the 2DEG system there are additional coherences due to 
the ground state 2DEG excitations, described by the last term. In addition to 
interactions among the photoexcited carriers similar to the undoped system 
[20], the latter term describes the scattering and correlations between MPs 
and spin- J, carriers, which lead to the relaxation of the photoexcited carriers 
due to MP emission and absorption. 

As can be seen from Eq. (40), the intraband scattering processes are described 
by density matrices of the form (X^X), (p%pf), {p^p'l), and (p\p\)- In the case 
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of the spin-f polarized 2DEG, the MP contributions are described by the 
density matrices (pf/3|) and (p|/5|), which vanish in the undoped system in 
the case of right-circularly polarized light and are analogous to the phonon- 
assisted intraband density matrices [20]. They relate an initial state consisting 
of a photoexcited electron or hole to a final state consisting of an electron or 
hole plus a MP and describe the effects of carrier scattering by MP emission 
or absorption. 

The number of independent density matrices can be reduced by noting that 
the e-h pair creation operators may be expressed in terms of exciton operators: 

pt U __LvX f e -i<l4k+k')l 2 /2, / 41 X 

e knl n -k'ml — ^njy 2-^^cpim e Oq y ,k-k^ V^-) 

obtained from Eq. (8). Similar to the DCTS, further reductions can be ob- 
tained by noting that, as discussed below, only many-body states with one 
valence band hole contribute to the above intraband density matrices in the 
case of the third-order non-linear optical response, and thus the density ma- 
trix (WWhh) can be neglected to this order. Furthermore, in the case of the 
spin-f polarized 2DEG excited with right-circularly polarized light, spin-j 
carriers are only created via the photexcitation, and thus the density matrix 
(ejeje|ej) contributes to higher order, similar to the undoped system. In this 
case, only states with one spin-j electron or hole contribute, and we obtain 
by using Eq. (41) and denoting by n e ± and n h ^ the number operators of the 
spin-j carriers that 

P%nn'i — P%nn'i n hl = ^2 e^ Q><q ^ ''^q'nm^q'-qn'm (42) 

and 

tfU'i = ftnn'inei = 4^E e^ (qxq,) ^ 2/2 ^, mri X q ,_ qmn , (43) 
Viv q , m 

The above expressions can be used to show that, in the case of spin-polarized 
2DEG (e.g. for filling factors v — 1/m (m = integer) or integer v) and right- 
circular polarization, the independent intraband density matrices have the 
form (X'X) and (p+X'X). More details will be presented elsewhere. For other 
filling factors, all density matrices that enter in Eq. (40) must be calculated 
to obtain the third-order response, with the exception of (wwhh) which con- 
tributes to higher order. 

To conclude this section, we note from the above equations of motion that the 
effects of the interactions on the non-linear optical response are described by 
the interband density matrix (Y) and the intraband density matrix (M n ). Due 
to the many-body nature of this strongly correlated system, approximations 
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are needed in order to calculate these interaction-induced effects. In the un- 
doped system, the DCTS cumulant expansions separate the coherent from the 
incoherent and the correlated from the uncorrelated contributions. In the case 
of the 2DEG, we introduced in [7] an analogous decoupling scheme based on 
a decomposition of the many-body semiconductor state \ip) that evolves from 
the ground state \G) according to the Schrodinger equation for the Hamilto- 
nian H tot (t). This method is summarized in the following section and then used 
to derive a factorization scheme for the density matrix equations of motion. 



5 Interaction Effects 

In this section we discuss a decomposition of the photoexcited many-body 
wavefunction into correlated and uncorrelated contributions, which we use in 
the next section to devise approximations for treating the interaction-induced 
density matrices (Y n ) and ([F n ,X^]). Similar to the DCTS, we expand in 
terms of the optical field and calculate the third-order polarization, which 
is expected to describe the non-linear optical signal when the photoexcited 
carrier density is small and the X-cold 2DEG correlations prevail. 

As in the theoretical approaches of [24,28], we note the one to one corre- 
spondence between the photon absorption/emission processes and the e-h 
pair creation/destruction. However, since here a 2DEG is present prior to the 
photoexcitation, when following the effects of the applied fields we count the 
number of valence band holes in a given state. Therefore, we use the shorthand 
notation 0-h, 1-h, 2-h ... to label the states, and it is clear that states with 
three or more holes do not contribute to the third-order non-linear polariza- 
tion [33]. We can then decompose the optically-excited state \ip) according 
to 

m = l^o) + l^i) + |V>2>, (44) 

where \ipi), % = 0, 1, 2, describes the contribution of the i-h states, with initial 
condition oo)) = 5i t o\G). From now on, we refer to operators that change 
the number of holes as interband and to operators that leave the number of 
holes unchanged as intraband. 



5.1 Linear response 

To lowest order in the optical field, only the 1-h state [^i) is photoexcited. 
We separate out the magnetoexciton contribution to this state by introducing 
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the decomposition 



IW = EiT-W + fc). ( 45 ) 

where \ip\L) is the {l-/i/2DEG*} contribution, defined by the condition (X^i) 
= 0, that describes the incoherent contributions due to the X-2DEG interac- 
tions. From now on, 2DEG* denotes excited 2DEG configurations. The exciton 
amplitude 

Pt = PQIVil) (46) 

coincides with the linear polarization, whose equation of motion is obtained by 
linearizing Eq. (31). As shown in [7], we obtain the equation of motion for |t/>il) 
by substituting the decomposition Eq. (45) into the linearized Schrodinger 
equation for and using the equation of motion for the linear polarization 

Pt- 

i 1 v i 

where the amplitude 

Pt = (m iL ) = (Yi\ip 1L ) (48) 

coincides with (Yj) to first order in the optical field. If we restrict to the LLO 
and LL1 states, Eq. (28) applies and we have that 

Pt = -Pt = P L . (49) 

The amplitude P t L describes the time evolution of the X+MP states that 
contribute to \Y) and corresponds to a X+MP coherence, analogous to the 
X+phonon coherence in undoped semiconductors [20,24]. 

To obtain the equation of motion for P L , we must consider the action of the 
Hamiltonian on the state \Y) introduced in section 3. Similar to Eq. (21) that 
defines the latter state, we introduce a new state \Z) orthogonal to all the X 
states as well as to \Y) as follows [7]: 

H\Y) = n\Y) + w (1*L - + \Z), (50) 

\ i — V\ i — Vq j 

where we obtain after some algebra [7] from the above orthogonality require- 
ments and after using Eqs. (10), (21), and (28) 

- (Y\H\Y) 

n= [ -^^,W=(Y\Y). (51) 
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By projecting the state (Y\ to Eq. (47) and using Eqs. (50), (28), and the 
orthogonality among the above states, we obtain the equation of motion for 

pL. 

id t P L = (ft - %i)P L + W (-^— - t^—) + Z L (52) 

V 1 - v \ l-voj 

where Z L = (Z\ip 1L ) and 7 is the dephasing rate. 

By continuing the above orthogonalization procedure, we create a Lanczos 
basis [29,42] of strongly correlated orthogonal states \Z^), where \Z^>) = \Y) 
and \Z^) = \Z), from the recursive relation obtained by acting with the 
Hamiltonian H on the previous state, and then orthogonalizing the result 
with respect to all the existing basis states [7,42]: 

H\Z {n) ) = Cl (n) \Z {n) ) + W {n) \Z {n ' 1] ) + \Z {n+l) ) (53) 

where 



Q(n) _ (Zn\H\Z n ) ( n ) _ (Z n \Z n ) 
(Z n \Z n ) \Z n -i\Z n -ii 



(54) 



By projecting the state (Z n \ to Eq. (47) we obtain the equations of motion for 
the coherences Z^ = (Z n \ipi L ): 

id t zt ] = (^ (n) - iln) Zt ] + W^Z^ + Zt +1) . (55) 

As discussed in [7], by taking the Fourier transform of the above system of 
equations, we can obtain a continued fraction expansion of the polarization, 
which describes a non-Markovian dephasing. This approach is analogous to 
the numerical calculations of the 2DEG dynamical structure factor [1,2]. The 
above hierarchy can be truncated when convergence is reached, which becomes 
more rapid with increasing damping rates or, in the case of an iV-electron 
system, after performing iV iterations. In the QHE literature, numerical cal- 
culations of the iV-electron spectral functions have been shown to extrapolate 
to the iV — > 00 limit for a relatively small iV [1,44]. Compared to an expan- 
sion in terms of noninteracting X+MP states, as in Eq. (26), the correlated 
Lanczos states are advantageous when the different momentum contributions 
are strongly coupled, or for obtaining a simple solution such as the average 
polarization model discussed below. 

In Fig.l we show the linear absorption spectrum obtained by retaining the 
states \X n ),n = 0, 1, and \Y). The effect of the higher Lanczos states, which 
describe the "bath" that leads to exciton dephasing and do not couple directly 
to the exciton, is taken into account by introducing the dephasing rate 7 of 
P L that describes the coupling between the exciton and the "bath" . This de- 
phasing is due to X+MP scattering (see Eq.(26)) and the emission of other 
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Fig. 1. Linear absorption spectrum with parameters as in [7,8,14] chosen to fit the 
experimental results at B = 7T (\/W = 3meV, T = F 1 = O.UmeV, 7 = 3.8meV, 
u = 0.34, i/i = 0). 

2DEG excitations. Noting the analogy between a MP and an X discussed in 
Section (3), one can make an analogy between the X-MP scattering described 
by Eq.(26) and the XX scattering in undoped semiconductors. As shown in 
[35] , in the case of magnetoexcitons the latter can be described by a dephasing 
rate for strong X-X interactions and leads to the average polarization model 
[16,30,35]. The main feature in Fig.l is the strong LL1 broadening and asym- 
metric lineshape, which is due to P L and cannot be obtained by introducing 
a polarization dephasing time. Even though the time evolution of the X+MP 
states, described by P L , determines the lineshape of the LL1 peak, it only 
plays a small role at the LL0 frequency. To interpret this difference we note 
that, as can be seen by using Eq. (26) to calculate the state \Y), the main 
contribution to P L comes from {1-MP + l-LL0-e + 1-LLl-h} four-particle 
excitations, while all other contributions are non-resonant with the LL0 and 
LL1 excitons. Even though P L couples equally to both X amplitudes Pq and 
Pi, it dominates the dephasing of P± since the above four-particle states have 
energy comparable to that of X\ in the case of an inter-LL MP. In contrast, 
Xq has significantly smaller energy, and thus the inter-LL MP plays a minor 
role in the broadening of the LL0 exciton peak. 

5.2 Second order processes 

In this section we consider the two-photon non-linear optical processes that 
lead to the photoexcitation of the 2-h state \ip 2 ) and the 0-h state |^ )- By sep- 
arating out the contribution of the states = Xjxj,\G), which describe 
a pair of non-interacting Xs, and Xjlipn), which describes a noninteracting 
pair of X and X+MP 1-h states, we arrive at the decomposition to 0(E A ) 

1 pL pL pL 
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where \ip2) describes the correlated XX and X-X+MP contributions and 
satisfies the equation of motion [7] 



(57) 

We note that, in the above equation of motion, there are no terms proportional 
to d(t) and therefore the decomposition Eq. (56) eliminates all contributions 
to \ip 2 ) that are proportional to the excitonic amplitudes ( whose time 
derivative is proportional to d(t)). 

In addition to the photoexcitation of the above 2-h many-body state, the 
two-photon processes of excitation of a 1-h state and then de-excitation of 
an e-h pair, possibly accompanied by the scattering of 2DEG excitations, 
leads to a second order contribution to the 0-h state IV'o)- We split the latter 
state into the contribution of the ground state \G), with amplitude (G\ip), 
where the 2DEG is not excited during the excitation-deexcitation process, and 
the photoexcited {0-/i/2DEG*} contribution. We further decompose the latter 
contribution into an uncorrelated part oc X^il), where the deexcited exciton 
Xi does not interact with the \4>il) carriers, and a correlated contribution |t/> ) : 

|^o) = (G\ip) \G) - ]T I±L Xl \^ L ) + |-0 O > + O(E'), (58) 
where the 2DEG* state |t/>o), (G'IV'o) — 0, satisfies the equation of motion [7] 

pL* pL pL*-y pL* -y- 



- E n 5vf z - E 7^7 ([*>*Jl - l G >- (59) 

U> l 1 _ ^JU _ v i') u> 1 ~ V V 

The first term in Eq. (59) describes the photo-excitation of the 2DEG via 
the second-order interaction-assisted process where the exciton X^, photo- 
excited with amplitude P^, scatters with the 2DEG into the state \Y V ), and 
then the exciton is deexcited with amplitude P t L . The above process leaves 
the system in a 2DEG* 0-h state. It is analogous to the photoexcitation of 
coherent phonons in undoped semiconductors, and dominates the inelastic 
light scattering spectra of the 2DEG. [45,47] The second term on the rhs of 
Eq. (59) describes the scattering of Xi with the carriers in j^) during its 
de-excitation. The rest of the terms describe the possibility to create 2DEG 
excitations by photoexciting an exciton whose hole then recombines with a 
2DEG electron (Raman process). In the case of the spin-f polarized 2DEG 
and right-circularly polarized light, the latter Raman process vanishes since 
there are no spin-j electrons in the ground state. 
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Below we use the above decompositions of the photoexcited many-body wave- 
function in order to describe the interaction-induced contributions to Eqs. 
(31) and (37). These decompositions provide a way of separating out the un- 
correlated from the correlated parts in a general correlated system, where a 
Hartree-Fock noninteracting state may not be an appropriate reference state 
as in undoped semiconductors. This separation also motivates a factorization 
of the corresponding density matrices that applies not only to undoped semi- 
conductors, but also to systems with strongly correlated ground states, where 
Wick's theorem does not apply. Similar to the DCTS, our method provides 
a systematic way of identifying the parts that can be factorized and the new 
intraband dynamical variables that cannot be expressed in terms of inter- 
band coherences due to the incoherent processes. Importantly, it allows us to 
treat both coherent and incoherent processes in strongly correlated systems, 
where the coupling between the photoexcited carriers and the "bath" (here 
the 2DEG) leads to new dynamics governed by slow/low energy "bath" col- 
lective excitations. Finally, our method raises the possibility of devising new 
approximations, obtained by projecting the many-body wavefunctions in an 
appropriate basis of strongly correlated states or operators, and then using this 
expansion to evaluate the correlated contributions to the density matrices. 

5.3 Intraband density matrix 

In this section we turn to the equation of motion for the density matrix (M), 
where M is any intraband operator that does not change the number of holes. 
Furthermore, we assume that (G\M\G) = 0, as is the case for the operator 
[Yj,Xj] discussed above. Substituting the above decomposition of the many- 
body state \if)) we obtain for the average value {ip\M\ip) by keeping terms up 
to second order in the optical field 

+ £T^<G|[*,M]fc) +T,^-(i;L\lM.x!]\G) +0(E% (60) 

where (M) c is the correlated contribution, given by 

(M) c = {G\M$ ) + $ \M\G) + fc|Mfc>. (61) 

This above result corresponds to an intraband density matrix decomposition 
into a factorizable part and a correlated part (M) c . The second term on the 
rhs is the coherent contribution, which similar to the undoped system can be 
expressed as a product of exciton polarizations, while the rest of the terms 
describe the incoherent contributions and 2DEG photoexcitation processes. 
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The above decomposition corresponds to a projection of the exciton states 

We can obtain the equation of motion for (M) c either from Eq. (61) by using 
the equations of motion Eqs. (47) and (59) and the orthogonality (G\M\G) = 
or from Eqs. (30) for (M) and Eqs. (31) and (47): 

id t {M) c - ([M, H]) c + i lM {M) c = £ P r ^' 



x [i(I> + I> - 7M )(X i ,|M|X ? -,) + (Yi,\X},M\G) - (G\MX lf \Y r 
+ E ' (r \ +7 ~ 7M) [Pe ($il\[M, X\,]\G) + (G\ [X,, M] fc> 

pL pL* pL pL* 

i'f l 1 - ^'H 1 - Vj') i'j> I 1 ~~ - V 3') 

+ E r^- I *?] 1 °) + E t^— < G I ' ^ fe> 
v i — j/ i — 

+ E t^— M]fc) + £ T^-^iLltX/,, m]\G) 
v l — ^i' j/ l — ^i' 

+ E Xt]M|G) - £ d ^t(G\M[X^ X),\\G). (62) 



The above equation of motion applies for any intraband operator M, (G\M\G) 
= 0, such as the operators that contribute to Eq. (40), and describes contri- 
butions due to the possible photoexcitation of a 2DEG coherence associated 
with the 2DEG* state M\G) as well as due to deviations from the factoriza- 
tion Eq. (60) induced by incoherent processes that involve the photoexcited 
carriers. The above equation of motion can be evaluated e.g. by expanding 
the operator M in the strongly correlated basis Z^)(Z^\, \X n )(Z^\, and 
\X n ){X m \ and a similar basis of 0-h 2DEG states discussed below. 



5.4 Interband density matrix 



We now consider the equation of motion for the density matrix (Y), where 
Y is any interband operator that creates an e-h pair with the simultaneous 
scattering of any number of other electrons or holes (e.g. the operator Yi nt , 
Eq. (26)). By using the above decomposition of the photoexcited many-body 
state we obtain after some algebra by keeping terms up to third order in 
the optical field the following decomposition of the interband density matrix 
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(Y) into correlated and uncorrected contributions [7]: 



I 1 % 1 

+ 1 E n P \f/ A ^l\[[y,xJ],xI]\g) + (y> c (63) 

2 ji# (1 - i/i){L - v v ) 

where 

(f) c = ^\G) P L + (Y>ijvl> + @o\Y\$il) + fc|^|^2>, (64) 
and we introduced the non-linear 1-h state 

\Anl) = ~ IVW + E T^^l^) - E t^-X^o). (65) 

The first term on the rhs of Eq. (63) describes the coherent X-X interaction 
and correlation effects analogous to the undoped system, discussed in the 
following section. The second term describes the contribution due to scattering 
of the optical polarization with the correlated intraband contributions, i.e. the 
intraband coherences and incoherent populations discussed above. The above 
two terms treat the effects of coherent X-X interactions and the intraband 
excitation and incoherent population processes. The third term in Eq. (63) 
describes X-X interactions accompanied by the shake-up of 2DEG excitations, 
while the last term describes the correlated contribution to the non-linear 
polarization. As can be seen from Eq. (64), the dynamics of the first two terms 
of the latter contribution is governed by the dephasing of the X+MP state 
\Y) while the last two non-linear terms are determined by the {l-/i/2DEG*} 
state \4>il) an d describe incoherent correlated interband contributions. The 
equation of motion for (Y) c can be obtained as above and describes the X 
dephasing. 



6 Coherent X X correlations 



In this section we make the connection between the above result for the 
coherent X-X interaction contribution to the non-linear polarization (first 
term on the rhs of Eq. (63)), described by the amplitude of the 2-h state 
{G\[Xi, Y]\ip 2 ) = ([Xi, Y]) +0(E 5 ), and the familiar expressions that describe 
X-X correlations in undoped semiconductors [23,24]. Using Eq. (26) and re- 
stricting to the first two LLs we obtain for B = [X 1: Y\ 



B 



q L 2 



oi(-q)A > _ qO i-0io(-q)X_ q io 0oi(q)^ q oi-0io(q)^ q io (66) 
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and [X , Y] = —B. As can be seen from the above equations, the coherent 
XX interactions are determined by 2-X density matrices of the form (XiXj), 
similar to the undoped system, where i,j describe excitons with finite total 
momentum. Using the decomposition Eq. (56) we obtain for any 2-h state \B) 
(e.g. the states \Y\ Xj]\G) that determine the coherent X-X contribution to 
the non-linear polarization) 



1 (R\ X X \ P L P L P L 



where B c = (_B|^ 2 ) describes the X-X and X-X+MP correlations. The first 
term in the above equation describes the familiar Hartree-Fock X-X inter- 
actions, while the second term describes the analogous interactions between 
the exciton X { and the {l-/i/2DEG*} state fc). By using Eq. (66) for B 
to calculate the overlap (G\B\XiXj) and noting that, in the X-phonon sys- 
tem, (G\BX}\iJjil) = 0, we reproduce the results obtained in [23,24] for the 
undoped system by using the cumulants. 

We now turn to the equation of motion for the correlated X-X amplitude B c . 
By projecting the state (B\ to Eq. (57), restricting to the LLO and LL1 states, 
and using Eq. (28) we obtain that 



L 





V\ 1 - Vq 



{B\X\\i, 1L ) {B\Xt$ 1L ) 



pL 

(68) 



By using the Lanczos recursive method [42] one can generate a basis of strongly 
correlated 2-h states similar to Section 5.1. By describing the effects of the 
higher Lanczos states by introducing a dephasing rate we recover the average 
polarization model results used to describe X-X correlations and biexciton 
effects in undoped semiconductors [16,30,35]. Thus the Lanczos method can 
be used to derive such a model, which, in the case of 2D magnetoexcitons, 
was shown in [35] to be a good approximation in the case of attractive or 
strong repulsive X-X interactions with appropriate range. The Lanczos basis 
leads to a continued fraction expression of the X-X amplitude B c , similar to 
the linear response of the 2DEG [7] , and is advantageous in the case of strong 
coupling between the different X momentum states, e.g. due to an antibound 
continuum resonance or a bound biexciton [35]. 
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7 Magnetoplasmon Dynamics in the coherent regime 

In this section we derive from the above theoretical formulation a generalized 
average polarization model that we then use to calculate the time-integrated 
three-pulse FWM signal in the coherent regime. In the latter regime, we ne- 
glect all incoherent contributions due to the photexcited carriers, which are 
described by \4>il)- As can be seen from Eq. (61), the intraband contribution to 
the FWM signal then comes from the coherent MP amplitude (G\[Y , Xj}]^) . 
Noting the relation 

$o\M n \G) =M* n + (S nl - 5 n0 )(i> \HAu n \G), (69) 

obtained by restricting to the first two LLs and using Eq. (34) and the eigen- 
value equation H\G) = 0, we obtain from Eq. (61) 

(M n ) c = M n + M* n + (5 nl - 5 n0 )(^ \HA0 n \G). (70) 

Similar to the X-X amplitude (B) c , the amplitude A4 n (t) describes the time 
evolution of the 0-h 2DEG* state M\\G) = X n \Y) discussed below. Using 
Eqs. (38) and (34) we obtain that 

\M) = M\\G) = - Ml\G) + HA0 \G) 

+ V 10 (l - z/ )Az> 1 |G) - V 10 (l - i/i)Ai>o|G) (71) 

which gives the following relation between Mi and M. : 

M(t) = M 1 (t) = - M (t) + (G\AP H\^ ) 

+ V 01 (l - u ){G\AuSo) - V 01 (l - z/i)(G|Az> |^o). (72) 

In the case of the spin-polarized 2DEG (QHE ferromagnet), only spin-f states 
are populated in the ground state and therefore Av n \G) = (see Eq. (19)). 
This is also the case if we neglect any fluctuations in the filling factor of the 
spin-j LLs in the ground state. We then obtain by using Eq. (71) that 

(G\M M j \G) = -(G\M Ml\G) = (G\MiM\\G) = ~{G\M Q M\\G) = W M 

(73) 

and, from Eq. (72), we obtain the relation 

M(t) = M 1 (t) = -M (t). (74) 

An explicit expression for the state \M) can be obtained by acting on the 
ground state with the intraband operator Eq. (40). Noting that the valence 
band is full in the ground state, the states that contribute to \M) are the MP 
states pl nm JG) and the states p e - mi p\ nma \G) and pt^P^JG) that vanish 
in the case of a spin-polarized ground state. As already known from inelastic 
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light scattering experiments, the MP states p\ nma \G) can be photoexcited 
in the real system due to the residual disorder, inhomogeneity, and valence 
band mixing [45,46,47,48]. For exaqmple, MPs with momentum much larger 
than the photon momentum, close to the maximum of the MP dispersion or 
close to the magnetoroton dispersion minimum q ~ l/l, have been observed 
in such experiments [45,46,47,48]. For the filling factors and photoexcitation 
conditions of interest here, the LLO — > LL1 MP will dominate, due to its longer 
lifetime and resonant contribution. 

To evaluate the equation of motion for M.{t) we also need the state H\M), 
which can be obtained by using the Lanczos approach 

H\M) = Q M \M) + \M) , Q M = , (M\M) = 0. (75) 

Neglecting \M) corresponds to the single mode approximation [1,2,50]. The 
above correlated basis set leads to a a continued fraction expansion for M.(lo) 
similar to the linear polarization calculation of Section 5.1 [7], the correlated 
XX amplitude calculation outlined in Section 6 [35], and the 2DEG dynamical 
structure factor calculations disucssed e.g. in [1,2]. 

Neglecting all incoherent contributions and using the above relations we obtain 
from Eq. (59) or Eq. (62) the equation of motion 

/ pL* pL* \ / pL pL \ 

idtM = (Q M - iim)M + W M -± —2— — ^- - — °- , (76) 



where 7m describes the damping of the 2DEG excitation \M). Neglecting the 
incoherent non-linear contributions to Eq. (64) we also that 

id t (Y) c = (n - < 7 )(y) c + w f-^- - -^-) , (77) 

\ 1 — V\ 1 — Uq J 

which describes the dephasing of the non-linear polarization similar to the 
case of the linear absorption. Substituting the above results into the equation 
of motion Eq. (31) we obtain for the third order non-linear polarization in the 
coherent regime 

id t P n - (Q n - lT n )P n + V nn ,{l - V n )P n , = 2 ^ E ^ P n P n + _ Sno) x 

V XX - ^3.) + (M+M-) (j^- - + (Y), 

V 1 — V\ I - Uq J V 1 — V\ I - VqJ 

(78) 

where n,n' — 0, 1, n ^ n', The first term on the rhs of the above equa- 
tion describes the Pauli blocking effects due to the coherent exciton density 
P n Pn- The second term describes the mean field X-X interactions, where 
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Vxx = (B\XiXq) is the only finite contribution to Eq. (67). The above terms 
reproduce the results of [53] if we restrict to the first two LLs. The third 
term in Eq. (78) describes the effects of the photoexcitated MP coherence 
M, whose time evolution is determined by the MP dynamics, while the last 
term describes the non-Markovian dephasing of the X polarization, which is 
governed by the time evolution of the X+MP states Using the above 
equations, we calculate the time-integrated FWM signal 

/+oo 
rft|P (t)+Pi(t)| 2 (79) 
-oo 

in the k s = k x + k 2 — k 3 direction, assuming right-circularly polarized optical 
excitation of the form E(t) = e ikl -%(f) + e ik2 ' r ^(HAr 12 ) + e ik3 '%(t + Ar 13 ), 
where £ p (t) is the Gaussian optical pulse. In the next section we discuss the 
results of this simple model that describes only coherent effects. The role of 
the MP in the population relaxation and the role of MP-assisted inter-LL 
coherences among the photoexcited carriers will be discussed elsewhere. 



8 Three pulse time integrated four—wave— mixing signal 

In this section we present the results of our numerical and analytical calcula- 
tions based on the model of Section 7 and using parameters of Fig. 1, obtained 
by fitting the experimental linear absorption as discussed in [7,8]. Below we 
analyze the time dependence of the TI-FWM signal as a function of two time 
delays: Ari3, the delay between pulses 1 and 3 also accessible in two-pulse 
FWM experiments, and Ari2, the time delay between pulses 1 and 2. We find 
that AT12 can give direct information on the MP dynamics that is not acces- 
sible as function of Ar 13 , where the time-dependence is similar to the results 
of [7]. We study below the FWM signal that comes from the three different 
non-linearities in Eq. (78), i.e. the PSF, X-X interactions, and MP coherence. 
Furthermore, we show that the correlation effects can be controlled by vary- 
ing the central photoexcitation frequency from LL1 toward LLO. Important 
is the fact that, as already seen in the linear absorption spectra (Fig. 1) and 
discussed in [5,7,8], the LL1 magnetoexciton decays much faster than the LLO 
magnetoexciton due to the X+MP dynamics described by P. For simplicity, 
when we discuss the numerical results below, we refer to exciton dephasing 
rates r and F 1 of P and P 1 respectively, with F 1 ^> r . In the experiment 
[14], the LL1 dephasing was measured to be of the order of a few hundreds of 
femtoseconds, while the LLO dephasing time was of the order of a few picosec- 
onds. Fig. 2 shows the PSF contribution to the TI-FWM signal, obtained by 
only retaining the first term on the rhs of Eq. (78), whose time dependence 
is determined by the dephasing of the LLO and LL1 magnetoexcitons. The 
PSF signal is maximum at Ar 12 = Ar 13 = and peaked along the directions 
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Fig. 2. PSF contribution to the TI-FWM signal when exciting (a) LLO and LL1 
equally (b) mostly LL1 (arbitrary units), with parameters as in Fig. (1). 

Ar 13 > 0, Ati2 = 0, Ar 12 > 0, Ar 13 = and Ar 13 = Ar 12 < 0. Fig. 2 compares 
the PSF contribution for optical excitation centered at LL1 (Fig. 2(b)) and be- 
tween the two LL peaks (Fig. 2(a)). Fig. 2 demonstrates a strong dependence 
of the PSF temporal behavior on the central photoexcitation frequency, which 
is due to the much faster dephasing of Pi as compared to Pq and the fact that 
the PSF contribution to Eq. (78) is directly proportional to the optical pulse. 
When the excitation frequency is centered on LL1, the PSF signal is domi- 
nated by the LL1 density contribution, while the LLO density is suppressed. 
In this case, the PSF signal decays very fast, as ~ 2r 1 (Ari2 + Ar^) in the 
Ati2 > 0, Ari3 > region and as ~ 2Fi(Ari 3 — 2Ari 2 ) in the Ati 3 > Atu, 
Ati2 < region. On the other hand, when both LLs are equally excited, the 
long-lived LLO exciton contribution dominates over the short-lived LL1 exci- 
ton contribution for sufficiently long time delays, and thus the signal decays 
more slowly as compared to the case of LL1 photoexcitation since To Ti. 
Moreover, beatings occur between the LLO and LL1 contributions, with fre- 
quency ~ Qi — fl , the difference in the energies between the two X peaks, 
and decay rate ~ r + We conclude that the PSF contribution to the 
TI-FWM experimental signal is characterized by a strong dependence on the 
central photoexcitation frequency, due to its proportionality to the optical 
pulse and the large difference in the dephasing of the two LL excitons due to 
the X+MP states. 

Fig. 3 shows the FWM signal due to the X-X interactions alone, obtained by 
only retaining the second term on the rhs of Eq. (78). The overall temporal 
profile is similar to the PSF signal, with strong oscillations as function of AT12 
with frequency — Q . By shifting the central photoexcitation frequency 
towards LLO, the XX signal along the Ari 3 > 0, Ari 2 = direction and the 
corresponding beatings are enhanced. 

The coherent MP contribution, obtained from the third term on the rhs of 
Eq. (78), is shown in Fig. 4. The new feature that distinguishes the MP con- 
tribution from the rest in the case of a long MP dephasing time is the strong 
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(a) 



(b) 




Fig. 3. TI-FWM signal due to the X-X interactions alone when exciting (a) LLO 
and LL1 equally, and (b) mostly LL1 (arbitrary units), with parameters as in Fig. 
(!)• 

signal in the directions At 12 < 0, Ati 3 = and At 12 = Ari 3 > 0, which is 
due to the photoexcitation of a MP with energy VLm comparable to the dif- 
ference between the LL1 and LLO exciton energies. This resonance between 
the MP energy and exciton splitting enhances the MP-induced signal, and 
leads to an optically-induced time-dependent hybridization of the LLO and 
LL1 magnetoexcitons, which is governed by the MP dynamics described by 
Ai and comes from the absorption of a photoexcited MP that couples the LL1 
to the LLO exciton as described by PiAi*. The above signal is mainly due to 
the resonant, oc P±Pq, contribution to the equation of motion Eq. (76), which 
describes a photoexcited inter-LL coherence that subsequently scatters into 
a MP 2DEG coherence via the X-2DEG interactions. Moreover, Fig. 4 shows 
strong oscillations as function of both time delays, which are mainly due to 
the non- resonant contribution to At, oc PoPq in Eq. (76), and are thus of 
quantum kinetic origin. By shifting the photoexcitation frequency from LL1 
toward LLO, P increases and thus the amplitude of these beatings, as well 
as the overall strength of the MP FWM signal is enhanced (see Fig. 4(a)). 
Of particular interest for studying the MP dynamics is the strong signal and 
oscillations in the directions Ar 12 < Ar 13 = and AT12 = Ar 13 > 0. The 
FWM signal in the above directions decays overall as ~ ^Jm, the MP de- 
phasing rate, and ~ 4r and displays oscillations with frequency ~ Qm, the 
MP energy, that decay with a rate of ~ 2r + Therefore, the above sig- 
nal can be used to extract from a three-pulse TI-FWM experiment the MP 
dynamics and dephasing rate. In contrast, the rest of the signal, along direc- 
tions Ar 12 > 0, At 13 = or Ar 13 > 0, Ar 12 = or At 12 = Ar 13 < decays 
overall as ~ 2r and shows oscillations with frequency ~ Qi — Q that decay 
with a rate ~ r + IY Finally, along the direction Ati 3 < 0, Ar 12 = 0, the 
signal decays overall with a rate ~ 4r and shows oscillations with frequency 
~ Qi — Q that decay with a rate ~ 3r + IY 

The main physical mechanism that gives the signal in the Ar 12 < 0, Ar 13 = 
direction is shown schematically in Fig. 5 and may be summarized as follows. 
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Fig. 4. Coherent magnetoplasmon contribution to the TI-FWM signal when excit- 
ing (a) LLO and LL1 equally, and (b) mostly LL1 (arbitrary units), with parame- 
ters as in Fig. (1). The MP correlations create a signal in the Ati 2 = AT13 > and 
Ati2 < = Ari3 directions that decays with the MP dephasing time = O.lmeV. 




Fig. 5. Non-Linear process that gives the MP FWM signal. This corresponds to an 
inelastic Raman process, assisted by the interactions between the photoexcited Xs 
and the 2DEG. 

The k 3 and ki pulses arrive at the same time and photoexcite an inter-LL 
coherence that subsequently scatters into a MP polarization due to the X- 
2DEG interactions. The k 2 pulse comes in later and creates an X polarization 
that scatters off the above decaying MP and gives a signal in the k x + k 2 — k 3 
direction. In particular, the MP scatters with the photoexcited X into an X 
state that then recombines to return the system in the ground state. It is 
interesting to note the similarity of the above process and the familiar one of 
coherent antiStokes Raman scattering [54] that involves phonons. Note here 
that, in the case of a spin-polarized 2DEG, the Raman process of excitation 
and de-excitation of the system with right-circularly polarized light is negligi- 
ble due the absence of spin-j electrons in the ground state. A MP excitation 
can however be excited due to the scattering of the photoexcited spin-]. X 
with the spin-f 2DEG. Obviously, if the MP dephases before the k 2 pulse ar- 
rives, there will be no FWM signal for this particular sequence of time delays, 
and therefore the above direction can be used to extract the MP dephasing 
dynamics from the experimental results. Furthermore, as discussed above, the 
~ Qm oscillations for Ar 12 = Ar 13 > or Ar 12 < 0,Ari 3 = directions 
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Fig. 6. TI-FWM signal with PSF + MP contributions when exciting (a) LLO and 
LL1 equally (b) mostly LL1 (arbitrary units) The parameters are as in Fig. (1), 
W M = W = V X x = 3 meV, j M = 0.1 meV. 

are of quantum kinetic origin and their decay gives important information on 
the MP coherent dynamics during time scales comparable to the MP period. 
Fig. 6 shows the TI signal when both the MP and PSF contributions are in- 
cluded. This figure shows a rich time-dependence, which is dominated by the 
MP contribution and can be used as a guide to extract from a three-pulse 
FWM experiment new information about the dephasing dynamics of the X 1 
and Xq excitons, as well as the MP dynamics during time scales comparable 
to the MP period. A comparison between the theory and such an experiment 
[14] as well as the role of the incoherent effects neglected in the above model 
will be presented elsewhere. 



9 Conclusions 



In summary, we discussed a recent theory [7,9,10] that provides a unified de- 
scription of the ultrafast non-linear optical response of magnetoexcitons in 
both doped and undoped semiconductors, including systems with a strongly 
correlated many-electron ground state, such as the 2DEG in the QHE regime. 
We discussed a method for describing the interaction contributions to the 
density matrix equation of motion, which gives the third-order non-linear po- 
larization measured in transient wave mixing and pump-probe experiments 
and established the connection with the DCTS treatment of X-phonon and 
X-X interactions in the undoped system [20,23,24]. Using a decomposition of 
the photoexcited many-body wavefunction into correlated/uncorrelated and 
excitonic/incoherent contributions, we obtained a factorization scheme for the 
density matrix that separates out the correlated from the uncorrelated con- 
tributions. Importantly, our method applies to systems with a strongly corre- 
lated ground state, populated by a 2DEG here, where previous factorization 
schemes, based on a Hartree-Fock ground state, need to be extended. Our for- 
mulation may be used to study the role of exciton-2DEG and exciton-exciton 
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correlations on the ultrafast non-linear optical response and the interplay be- 
tween coherent and incoherent effects. Our expansion in terms of the optical 
field is valid for sufficiently short pulses and/or weak excitation conditions, 
where the correlations are most pronounced. We describe the role of the long- 
lived collective excitations of a strongly correlated cold electron gas, which is 
present prior to the optical excitation. We also presented a number of results 
that describe the Pauli exchange and interaction effects between excitons and 
magnetoplasmon collective excitations. 

We applied our theory to the case of the photoexcitation of the 2DEG with 
three right-circularly polarized optical pulses. Our numerical solution for the 
time-integrated three-pulse FWM signal demonstrates a signature of the col- 
lective 2DEG excitation dynamics that, for long MP dephasing time, can be 
distinguished from the Pauli blocking and exciton-exciton interaction contri- 
butions. In the case of interest here, the relevant 2DEG collective excitations 
are the long-lived inter-LL magnetoplasmons, which dress the photoexcited 
magnetoexcitons and lead to polaronic-like effects and strong non-Markovian 
dephasing and quantum kinetic effects. We showed that such effects domi- 
nate the time delay and frequency dependence of the transient FWM signal in 
the coherent regime for time scales comparable to the inverse MP frequency. 
FWM spectroscopy using femtosecond optical pulses provides both the time 
and the frequency resolution necessary to access a new regime of 2DEG and 
QHE physics. Our theory allows us to study the experimental signatures of the 
2DEG quantum dynamics. We predicted, in particular, strong non-Markovian 
effects, leading to an asymmetric LL1 exciton lineshape and to strong oscil- 
lations as function of time delay of quantum kinetic origin that are governed 
by the MP dynamics. We also found an optically-induced hybridization of the 
LLO and LL1 excitons, with the time-dependent coupling determined by the 
MP dynamics. A comparison of our results with recent experiments [14] is 
currently in progress and will be presented elsewhere. 

The above correlation-induced non-Markovian dynamics can be controlled by 
tuning the central frequency of the optical excitation between the two lowest 
LLs, which changes the coherent admixture of the two MP-dressed magne- 
toexcitons. FWM experiments using a sequence of optical pulses with different 
polarizations and central frequencies provide new ways for accessing the very 
early dynamics of the strongly correlated 2DEG, during time scales shorter 
than the duration of the interaction processes. This opens up a new window 
into non-equilibrium Quantum Hall effect physics and collective effects that 
are just now starting to be explored. In particular, interesting regimes include 
(1) filling factors away from v — 1, where skyrmions [2] become important in 
one-sided modulation doped quantum wells and are expected to dominate the 
magnetoexciton dephasing, and (2) fractional filling factors, where intra-LL 
collective excitations will dress the magnetoexciton leading to polaronic effects 
[44]. 
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